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ABSTRACT 

SN 2006bt displays characteristics unlike those of any other known Type la supernova (SN la). 
We present optical light curves and spectra of SN 2006bt which demonstrate the peculiar nature of 
this object. SN 2006bt has broad, slowly declining light curves indicative of a hot, high-luminosity 
SN, but lacks a prominent second maximum in the i band as do low-luminosity SNe la. Its spectra 
are similar to those of low-luminosity SNe la, containing features that are only present in cool SN 
photospheres. Light-curve fitting methods suggest that SN 2006bt is reddened by a significant amount 
of dust; however, it occurred in the outskirts of its early-type host galaxy and has no strong Na D 
absorption in any of its spectra, suggesting a negligible amount of host-galaxy dust absorption. C II 
is possibly detected in our pre- maximum spectra, but at a much lower velocity than other elements. 
The progenitor was likely very old, being a member of the halo population of a galaxy that shows no 
signs of recent star formation. SNe la have been very successfully modeled as a one-parameter family, 
and this is fundamental to their use as cosmological distance indicators. SN 2006bt is a challenge to 
that picture, yet its relatively normal light curves allowed SN 2006bt to be included in cosmological 
analyses. We generate mock SN la datasets which indicate that contamination by similar objects 
will both increase the scatter of a SN la Hubble diagram and systematically bias measurements 
of cosmological parameters. However, spectra and rest-frame i-band light curves should provide a 
definitive way to identify and eliminate such objects. 

Subject headings: supernovae: general — supernovae: individual (SN 2006bt) 



1. INTRODUCTION 



A general model of Type la supernov a (SN la) explo- 
sions was established decade s ago (e.g., IW helan fc IbenI 

IColgate fc McKcc 1969,; 
a CO white dwarf accretes 



1971 



Arnetd 



Nomoto et al] 119841: 
1982f ). In tiiis model, 



material from a companion star until it reaches the den- 
sity for carbon ignition, causing a runaway nuclear reac- 
tion to propagate through the star. During this process, 
~O.6M0 of the star burns to ^^Ni which powers the lu- 
mino sity of the SN by decayin g to ^^Co and then to ^^Fe 
Ce.g.. [Colgate fc McKeel[l96l . 

This scenario predicts a close correspondence be- 
tween the amount of ^^Ni, the peak luminosity of 
the SN, and the temperature of the ej ecta wh i ch de - 
termines the ionization of the ejecta (jArnettI Il982f ). 
These predictions have been verified by observations of 
normal SN e la as well as both high-luminosity (e.g. , 
SN 1991T: iFifopenko et al.lll992bl : IPhpi ps et al."l992l) 
and low-luminosity (e.g., SN 1 991bg; [Filippcnko et al] 
Il992at iLeibundgut et all Il993f l SNe la. ' These rela- 
tionships pred ict a spectral sequence from cool to hot 
photospheres ijHachinger et al.l 120081) which have als o 
been empirically determined (e.g.. iNug ent et al.lfl99l . 
The width- luminosity relationship (WLR) between light- 
curve shape and luminosity that empowers cosmologists 
to use SNe la a s accurate and precise distance indicators 
()Phillipslll993f ) is directly related to the ionization (and 

^ Harvard-Smithsonian Center for Astrophysics, 60 Garden 
Street, Cambridge, MA 02138. 

^ Clay Fellow. Electronic address rfoley@cfa.harvard.edu . 

^ Department of Physics, Harvard University, 17 Oxford Street, 
Cambridge, MA 02138. 

^ Department of Astronomy, University of California, Berkeley, 
CA 94720-3411. 



therefore the temperature) of the SN photosphere, specif- 
ically when doubly ionized F e-group elements recombine 
to singly ionized ions IKasen fc Wooslevl[200l . 

Except for a few e xamples of ver y pecuHar SNe la 
such as SN 2000cx jLi et al.' '200ll ), the SN 2002cx- 
like class of objec ts (jLTet al. 2003; F'hillips et al]|2007t 
Foley et a^ll2009a^■ hv d rogen -rich objects (jHamuv et al.l 
2003t lAldering et all I2006D. and potential super- 



Chan drasekhar objects (|Howell et al.ll2006l : iHicken et al.l 
l2007t l. SNe la are observationally well described by their 
light-curve shape alone. There are methods of further 
subclassifying objects based on spectral ch aracteristics 
Benetti et all l2005t [Branch et al.l l2006f) . iFolev et all 



200a) an d lWang et al.l (|2n09D showed that incorporating 



spectral features into distance measurements from light- 
curve fitting can improve SN distance measurements, but 
the former study was based on a small sample and it 
is unclear if the improvement seen in the latter study 
is intrinsic to the explosions or related to the SN envi- 
ronment. Most physical parameters and observables of 
SNe la, such as ^^Ni mass, peak luminosity, light-curve 
shape, photospheric temperature, and ejecta ionization, 
can be parameterized by a single number (excluding clear 
outliers like the objects mentioned above); hence, know- 
ing any one of the parameters or observables will tell you 
the rest. 

SN 2006bt was discovered bv lLee fc Li' (2006) on 2006 
April 26.5 (UT dates will be used throughout this pa- 
per) at (R.A., Dec.) = (15:56:30.526, +20:02:45.34) 
(J2000), 44.4" west and 22.9 " south of CGCG 1 8-013 
the presumed host galaxy. iFilippenko fc FolevI (|2006l ) 
obtained a spectrum one night after discovery and de- 
termined that SN 2006bt was a SN la with some slight 
peculiarities. Photometric observations of the SN were 
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then performed using the 1.2 m telescope located at 
the Fred L. Whipple Observatory (FLWO), and light 
cu rves derived from th ese observations were published 
bv lHicken et all (l2009al) (hereafte r H09a). Fitting these 
hght curves. iHicken et all ()2009b( l (hereafter H09b) de- 
termined that SN 2006bt peaked in the B band on 
2006 May 1.4, with Ami5(B) = 1.09 ± 0.06 mag and 
Ami5(y) = 0.54 ± 0.04 mag, indicating that it had a 
slightly broad light curve and should thus have a sl ightly 
higher than average peak luminosity (jPhillipsH lQQS) . Fit- 
ting the light curves wit h more sophistic ated techniques, 
they found the SALT2 ([Gu7eral2?Ky^ and MLCS2k2 
(with Rv = 1.7: lJha et al.l[2007[ ) fit parameters of xl = 
0.069±0.104, c = 0.162±0.008 and A = -0.325±0.052, 
Ay = 0.428 ±0.053 mag, respectively. All light-curve fits 
produced reasonable values of ^nd there was no obvi- 
ous reason to suspect any peculiarity from quality cuts. 
Both SALT2 and MLCS2k2 indicate that SN 2006bt had 
a high luminosity and a red color (which MLCS2k2 inter- 
prets as a significant amount of host-galaxy extinction, 
but SALT2 does not attempt to distinguish between in- 
trinsic color differences and reddening from dust). Us- 
ing these light-curve fits, SN 2006bt was includ ed in a 
cosmological analysis of a large set of SN data ()H09b[ ). 
SN 2006bt was not an outlier in the Hubble diagram and 
its inclusion in the analysis did not significantly change 
the measured cosmological parameters. 

We present our data on SN 2006bt in Section [21 per- 
forming a detailed analysis of the photometry and spec- 
troscopy in Sections [3] and [H respectively. In Section El 
we investigate the SN environment. We discuss possible 
physical interpretations of SN 2006bt as well as implica- 
tions for SN la cosmology in Section [6l Our results are 
summarized in Sectional 

2. OBSERVATIONS AND DATA REDUCTION 

SN 2006bt was followed photometrically by the 1.2 m 
telescope located at FLWO and it s UBVri light curves 
were originally presented by lH09al . We reproduce those 
light curves in Figure [T] Because our conclusions are 
extremely sensitive to the nature of the light curves, we 
have re-reduced the data and compared our photometry 
of the field stars used for photometric calibration against 
photometry of the stars from SDSS in ri and V^; all 
measurements were found to be within 0.01 mag of the 
photometry published bv lH09a . 

We also obtained several low-resolu tion spectra of 
SN 2006bt with the FAST spectrograph (jFabricant et all 
[1998) on the FLWO 1.5 m telescope (as part of the 
CfA SN monitoring program; Blond in et al., in pre p.), 
the Kast double spectrograph (.Miller fc Stond [19931 ) on 
the Shane 3 m telescope at Lick Observator y, and the 
Low Resolution Imaging Spectrometer (LRIS: 10ke et al.l 
[1995) on the 10 m Keck I telescope; a journal of observa- 
tions can be found in Table [U Standard CCD processing 
and spectrum extraction were performed with IRAF^. 
The data were extracted using the optimal algorithm of 

^ http:/ /www. sdss.org/dr7/algorithms/ 
sdssUBVRITransform.html#Lupton2005 

^ IRAF: the Image Reduction and Analysis Facility is distributed 
by the National Optical Astronomy Observatory, which is operated 
by the Association of Universities for Research in Astronomy, Inc. 
(AURA) under cooperative agreement with the National Science 
Foundation (NSF). 
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Fig. 1.— UBVri hght curves of SN 2006bt shifted by 1.2, 0.27, 
0, —0.73, and —1.8 mag, respectively. Also plotted are the light 
curves of SNe 2005mz (Ami5(iJ) = 1.96 mag; solid lines) and 
2006ax (Ami5(B) = 1.08 mag; dashed lines). The light curves of 
SNe 2005mz and 2006ax have been shifted to match the peak of 
SN 2006bt in each band. The UBV light curves of SN 2006bt are 
very similar to those of SN 2006ax, while the r and i bands show 
significant differences near the "shoulder" in the r band and the 
"trough" and second maximum in the i band. [See the electronic 
edition of the Journal for a color version of this figure.] 

iHornd (|198(3t ). Low-order polynomial fits to calibration- 
lamp spectra were used to establish the wavelength scale, 
and small adjustments derived from night-sky lines in the 
object frames were applied. We employed our own IDL 
routines to flux calibrate the data and remove telluric 
lines using the well-e xposed continua of spectrophoto- 
metric standard stars (jWade fc Horndll988l : iFolev et all 
120031 ). Our spectra of SN 2006bt are presented in Fig- 
ure [1 

3. PHOTOMETRIC ANALYSIS 

3.1. Light Curves 

The lig ht cur ves of SN 2006bt have already bee n pub- 
hshe d by lH09al . The light curves were fit by both[H09i 
and [H09b|, finding Ami5(B) = 1.09 ± 0.06 mag and 
Ami^iV) = 0.54 ± 0.04 m ag. They similarly found 
the SALT2 (iGuv et al.| [2007l) and MLCS2k2 (with Ry = 
1.7);[JliaeLalJl20Q7D fit parameters ofxl = 0.069±0.104, 
c = 0.162 ± 0.008 and A = -0.325 ± 0.052, Ay = 
0.428 ±0.053 mag, respectively. The light-curve fits indi- 
cate that SN 2006bt has a slightly slow decline, and a red 
color (SALT2) or relatively high extinction (MLCS2k2). 
The MLCS2k2 fit has a reduced of -1.36 for 98 sep- 
arate observations, including pre-maximum data in mul- 
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TABLE 1 

Log of Spectral Observations 







Telescope / 


Exposure 




Phase^ 


UT Date 


Instrument 




Observer'' 


Q Q 
— O.O 


zUUd Apr. Z 1 .0 


ivecK l/ljillo 


oUU 




—3.0 


zUUd Apr. Z0.4 


Lick/Kast 


1800 


DW, JS 


— Z.KJ 


iSUUD Apr. zy.'i 


r ijVV w / ri\o 1 




VV D 


0.0 


2006 May 1.4 


FLWO/FAST 


1500 


WB 


0.9 


2006 May 2.3 


FLWO/FAST 


1500 


WB 


1.9 


2006 May 3.3 


FLWO/FAST 


1500 


TG 


2.9 


2006 May 4.3 


FLWO/FAST 


1800 


TG 


3.9 


2006 May 5.3 


FLWO/FAST 


1500 


WP 


4.0 


2006 May 5.4 


Lick/Kast 


1800 


JS, MM, RF 


5.9 


2006 May 7.3 


FLWO/FAST 


1500 


WP 


6.9 


2006 May 8.3 


FLWO/FAST 


1500 


WP 


17.9 


2006 May 19.3 


FLWO/FAST 


1800 


MC 


34.9 


2006 Jun. 5.4 


Lick/Kast 


2100 


JS, MM 


56.0 


2006 Jun. 26.4 


Lick/Kast 


2400 


DW, JS, TS 



^ Days since B maximum, 2006 May 1.4 (JD 2,453, 
AF = A. Filippcnko, DW = D. Wong, JS = J 
Calkins, MM = M. Moore, RF = R. Foley, TG = 
Steele, WB = W. Brown, WP = W. Peters. 



856.9). 
Silverman, MC = M. 
T. Groner, TS = T. 
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Fig. 2. — Optical spectra of SN 2006bt. The phase relative to 
B maximum is noted to the left of each spectrum. The zero flux 
level for each spectrum is indicate by the dashed line. 



tiple passbands. The SALT2 templates used by lH09al 
do not cover the late-time photometry and the best fit 
has a reduced of ~2 with 70 observations. However, 
estimates of the time of maximum and the best-fit light- 
curve models using either method are quite comparable. 
An analysis of the light-curve fits alone would lead one to 
conclude that SN 2006bt is not a particularly discrepant 
object and can be included in any cosmological analysis. 



However, the host galaxy, CGCG 108-013, is an SO/a 
galaxy (with d = l39Mpc) and SN 2006bt has a 
projected galactocentric distance (PGCD) of 33.7 kpc 
(equivalent to 7.8i?off, where i?off = 6.4" is defined by 
the Petrosian radius where 50% of the flux is contained 
within a circle of radius Ros] see Section [S]), making 
SN 2006bt have the fourth largest PGCD of SNe la in the 
CfA3 sample of 185 SNe la. The SN environment is thus 
expected to have little dust (excluding possible circum- 
stellar dust). Additionally, H09b shows the relationship 
between host morphology, PGCD, and light-curve fit pa- 
rameters. In the CfA3 sample, SN 2006bt has among the 
largest values for xl, c, and Ay, and the smallest values 
of A for all SNe hosted in early-type galaxies. Further- 
more, SN 2006bt is a clear outlier in the trend of smaller 
c or A y with increasing PGCD (see Figures 18 and 19 in 

MM- 

From the light curves, we can see that SN 2006bt 
evolves slightly differently from other SNe with simi- 
lar S-band decline rates. In Figure [1] the light curves 
of SN 2006bt are compared to those of SN 2006ax, an 
object with Ami5{B) = 1.08 ± 0.05 mag and minimal 
{Ay < 0.04 mag) host-galaxy extinction. These ob- 
jects have very similar UBV light curves; however, the r 
and i light curves have some slight differences. In the r 
band, SN 2006bt lacks the "shoulder" at 15-20 days after 
maximum brightness seen in average and high-luminosity 
SNe la. In the i band, SN 2006bt has a shoulder, but 
lacks the deep "trough" and subsequent prominent sec- 
ond maximum seen in most SNe la. 

In Figure [1] the light curves of SN 2006bt are also 
compared to those of SN 2005mz, a fast-declining SN la 
(Ami5(i3) = 1.96±0.14 mag) with moderate host-galaxy 
extinction {Ay = 0.27 ±0.09 mag). Although this object 
lacks a shoulder in either the r or i bands, it declines 
much faster in all bands. 

3.2. Color Curves 

Figure [3] presents the color curves of SNe 2005mz, 
2006ax, and 2006bt We have corrected the color curves 
of SN 2005mz for host-galaxy extinction as determined 
by MFCS fits {Ay = 0.27 mag and Ry = 1.7), and 
all objects for Milky Way extinction as determined by 
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Fig. 3.— B-y, y-r, and y-i color curves of SN 2006bt. The 
color curves of SNe 2005mz (Ar?ii5(_B) = 1.96 mag; corrected for 
Ay = 0.27 mag; solid lines) and 2006ax (Ami5(iJ) = 1.08 mag; 
dashed lines) are overplotted. We also plot the color curves of 
SN 2006ax reddened by Ay = 0.428 mag (with Ry = 1.7; dot- 
ted lines) to match the mea sured reddening determined by the 
MLCS2k2 fit to SN 2006bt IIHOQal) . The color curves are very 
different from those of both a normal SN la with a similar de- 
cline rate (SN 2006ax) and a low-luminosity SN la with a similar 
maximum-light spectrum (SN 2005mz). Furthermore, although the 
reddened B — V color curve of the normal SN la is similar to that 
of SN 2006bt, the V — r and V — i color curves are quite different 
in the ranges 10 < t < 25 days and < t < 30 days relative to 
maximum, respectively. It is worth noting that after the extinc- 
tion correction, SN 2006bt has a similar color to a normal SN la at 
t = 35 days after maximum brightness, when the color s of unred- 
dened SNe la are very similar regardless of decline rate IILira et al.l 
119981) . [See the electronic edition of the Journal for a color version 
of this figure.] 



iSchlegel et all (|1998f) . which indicates E{B - V) ^ 
0.051 mag for SN 2006bt. The curves of SN 2006bt gener- 
ally fall between the curves of SNe 2005mz and 2006ax for 
t < 35 days and above both curves for t > 35 days. Red- 
dening the color curves of SN 2006ax by Ay — 0.428 mag 
with Ry = 1.7, the best-fit value for SN 2006bt from 
lH09bl using MLCS, theB-V color curves of SNe 2006ax 
and 2006bt are very well matched. This i s consistent with 
the "Lira relation" (i Phillips et "anil999f) . which uses the 
late-time colors of SNe la to predict their host-galaxy 
reddening. However, this same correction does not ac- 
count for the differences in the V — r and V — i color 
curves of SNe 2006ax and 2006bt. Given the shapes of 
their color curves, reddening alone cannot account for 
the differences between these objects. Furthermore, the 
color curves of SNe 2005mz and 2006bt differ by more 
than reddening alone. 



Considering the PGCD, the early-type host, the non- 
standard color curves, and a lack of Na D absorption in 
any spectrum (see Section 2]), we believe the MLCS-fit 
extinction of Ay = 0.428 mag to be unrepresentative of 
the true host-galaxy extinction. Given these factors, the 
extinction is most likely very small, and for the remainder 
of our analysis, we assume that there is no host-galaxy 
extinction. 



3.3. Luminosity 

Since SN 2006bt has slightly odd light and color curves 
as well as a host-galaxy extinction very different from its 
MLCS fit extinction (which indicates that SN 2006bt is 
not well modeled by the training set of MLCS), it is 
worth investigating the possibility that SN 2006bt does 
not follow the WLR. The MLCS model would indicate 
that SN 2006bt had a peak absolute magnitude in V that 
is 0.428 mag (equivalent to the fit value for Ay) brighter 
than its true value. 

We plot the relationship between Ami5(i?) and My i n 
Figure H] for a subsample of the CfA3 sample (|H09al ). 
Using the extinction derived from MLCS, SN 2006bt 
is slightly overluminous relative to the relationship fit 
by iPhillips et al.l ()1999D . but it is within the scatter 
around the relationship. With no extinction correction, 
SN 2006bt is slightly underluminous for its light-curve 
shape; however, it is not a significant outlier. 

Despite an atypical set of light curves, SN 2006bt is 
not a significant outlier to the WLR using either no ex- 
tinction or t he ML CS-determined extinction value. For 
SN 2006bt, lH09bl measured distance moduli of fi — 
35.981, 35.944, 35.784, and 35.993 mag with distance fit- 
ters SALT, SALT2, MLCS with Ry = 3.1, and MLCS 
with Ry = 1.7, respectively. Using the CMB-corrected 
redshift for CGCG 108-013 of zcmb = 0.0330 and Hq = 
65 km s~^ Mpc~^ (to match the analysis of H09b), 
CGCG 108-013 has a Hubble-flow luminosity distance 
modulus of /i = 35.968 mag^. Therefore, three out of 
the four distance fitters provide values within 0.03 mag 
of the expected value if CGCG 108-013 is in the Hubble 
flow. We can also measure the distance modulus with 
the formula fj, = Fpoak - Myp^a^k - ^y.host - ^v mw- Us- 
ing the values measured and assumed bv IH09E we find 
H = 36.254 and 35.826 mag for Ay = and 0.428 mag, 
respectively. These values differ from the Hubble-flow 
distance modulus by 0.286 and —0.146 mag, respectively, 
corresponding to the offset from the WLR in Figure 21 
Obviously, fitting all light curves of SN 2006bt provides a 
better estimate of the distance to SN 2006bt than using 
only y-band data. 

Although SN 2006bt appears to have a somewhat odd 
light curve, it is still a relatively good standardizable can- 
dle. SN 2006bt is both intrinsically faint and red for its 
light-curve shape. All light-curve fitters correct for its 
red color by effectively brightening its apparent magni- 
tudes. Since intrinsic color variation and dust extinction 
are not orthogonal, deviations from t he model in MLCS 
are regarded as dust extinction (see iConlev et al.ll2007l 
for a discussion of how MLCS and SALT treat color). 

^ The host of SN 2006bt, CGCG 108-013, is a likely member of 
the Hercules supercluster and may have a large peculiar velocity. 
A peculiar- velocity uncertainty of 1000 kms~ corresponds to an 
uncertainty in the distance modulus of 0.012 mag. 
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Fig. 4. — Relationship between peak absolute magnitude {My) 
and decline rate (Ami5(i3)) for a subsample of the CfA3 dataset 
with Ay < 1 mag, E{B - V)mw < 0.5 mag, z > 0.015, and 
no major spectroscopic or photometric peculiarities (for e xample, 
SN 2002cx was excluded) . To be consistent with IH09bl we set 
Ho = 65 km s-l Mpc"!. The size of each point is inversely 
proportional to its uncertainties. A few points have uncertainties 
plotted to show the relationship. Overplotted is the quadrat i c rela- 
tionship between Amir^{B) and My given bv lPhillips et all II1999I ) 
(Mv(Ami5 = 1.1 mag) = —19.48 mag; solid lines) as well as that 
relationship shifted by an amount corresponding to the intrinsic 
scatter in the relationship (dotted lines). SN 2006bt is plotted 
both with its true value for My (blue point with solid error bars) 
and corrected by the MLCS fit value for Ay (red point with un- 
certainties indicated by the vertices of the surrounding diamond). 
Both points are consistent with the relationship. [See the electronic 
edition of the Journal for a color version of this figure.] 

We discuss the imphcations of including SN 2006bt-like 
objects in a cosmological analysis in Section 

4. SPECTRAL ANALYSIS 

In Figure [21 we show our collection of spectra for 
SN 2006bt. In this section, we examine the spectra in 
detail. 

4.1. Supernova Redshift 

SN 2006bt was located between CGCG 108-013 and 
2MASX J15562803-I-2002482, two objects at similar, but 
significantly different, recession velocities {cz = 9628 ±55 
and 13, 894±30 kms~^, respectively). To determine the 
most likely host galaxy, we cross-correlated all spectra 
obtained a tt<7 days with a libra ry of SN la spectra us- 
ing SNID (jBlondin fc Tonrvl [20071) . We determined two 
likely redshifts for each spectrum: the redshift from the 
best-fitting template and the average redshift from all 
SN 1991bg-like objects (since SN 2006bt has a very sim- 
ilar spectrum to this class of objects; see Section r4.2p . 
The best-fitting template spectrum was of SN 1986G 
for all but one spectrum. Performing these fits, we 
find C2SN = 9397 ± 109 and 8682 ± 73 kms'^ for the 
best-fit template and SN 1991bg-likc templates, respec- 
tively. Considering the velocity dispersion of the galaxy 
and the peculiarity of SN 2006bt, we believe that these 
values are consistent with SN 2006bt being associated 
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Fig. 5. — The pre- maximum-light spectrum of SN 2006bt com- 
pared to those of SNe 1986G {Amis{B) = 1.65 mag), 2006D 
{Ami5{B)approxl.35 mag) and 2006ax (Ami5(_B) = 1.08 mag) 
at a similar phase. The spectrum of SN 1986G has been dered- 
dened by E{B - V)mw = 0.115 mag and E{B — F)host = 0.4 mag. 
The spectrum of SN 2006bt shows many similarities to that of 
SN 1986G including strong Si II A5972 and the presence of Ti II. 



with CGCG 108-013 and inconsistent with SN 2006bt 
being associated with 2MASX J15562803+2002482. For 
the remainder of the paper, we assume the redshift of 
SN 2006bt is that of CGCG 108-013. 

4.2. Pre-Maximum Spectrum 

Our first spectrum was obtained 3 days before i3-band 
maximum. Using this spectrum, iFilippenko fc FolevI 
(|2006[ ) noted that SN 2006bt was somewhat peculiar. 
In Figure [5l we present this spectrum as well as spec- 
tra of SNe 1986G (Ami5(S) = 1.65 mag) and 2006ax 
(Arni5(i?) = 1.08 mag) at a similar phase. Despite 
having the same value of Ami5(i3), SNe 2006ax and 
2006bt exhibit very different spectra at this epoch, with 
SN 2006bt having a much stronger Si II A5972 fea- 
ture and a depression at ~ 4100 A corresponding to 
Ti II. These features are indicative of a relatively cool 
photosphere and are hallmarks of th e low-luminosity 
SN 1991bg-lik e class of SNe la (e.g., IFilippenko et aij 
[l992a; Leibun dgut et al.lll993D . Although there is some 
dispersion in t he spectra of obje c ts with the same light- 
curve shape ffBenetti et al.ll2004t iMatheson et al.l[200&) . 
differences of this magnitude have never before been seen 
within normal SNe la of the same light-curve shape. 

The earliest spectrum of SN 2006bt is more similar 
to the spectrum of SN 1986G, a low-luminosity SN la 
and a member of t he SN 1991bg-like class of objects 
(iPhillips et al.l[l98l . SN 1986G has stronger Si II A5972 
and Ti II features, indicating that SN 2006bt is slightly 
hotter than SN 1986G. (However, we do caution that the 
uncertain reddening correction of SN 1986G may cause 
the relative strengths of its features to be incorrect.) 

At this phase, the minimum of the Si II A6355 absorp- 
tion feature is blueshifted by 12,500 kms~^, which is the 
same as SN 1986G and within the typical range of slower 
declining SNe la (although it is slightly larger than the 
velocity of SN 2006ax at that phase). The minimum of 
the Ca II near-infrared triplet (seen in Figures [2l and [9]) 
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Fig. 6.— The t = +35 day spectrum of SN 2006bt com- 
pared to those of SNe 1986G {Ami5(S) = 1.65 mag) and 2006ax 
(Ami5(i3) = 1.08 mag) at a similar phase. The spectrum of 
SN 1986G has been dereddened by E(B - V)mw = 0.115 mag 
and E{B — y)host = 0.4 mag. At this phase, all three SNe have 
similar spectra. 

is blueshifted by 15,900 kms""'^, which is also normal. 
There is no obvious high-velocity structure associated 
with either feature. 

Despite having a slow optical decline rate, SN 2006bt 
has a pre-maximum spectrum most similar to those of 
fast decliners. SN 2006bt has a relatively cool photo- 
sphere, which may explain the lack of a shoulder in r and 
a prominent second maximum in i (Kascn 2006), but is 
contr adictory given the slow decline (|Kasen fc WooslevI 
[200l . 

There is no indication of Na D absorption at either 
zero velocity (c onsistent with the low value of E{B — 
V) measured bv lSchlegel et aLl fTgOS) or the host-galaxy 
recession velocity in any of our spectra, indicating a small 
or negligible extinction. 

4.3. Late-Time Spectrum 

In Figure [6l we present our t = 35 day spectrum com- 
pared to spectra of SNe 1986G and 2G06ax at a simi- 
lar phase. The differences between SN spectra become 
smaller with time, and by t = 35 days, the spectrum of 
SN 2006bt is very similar to spectra of both SNe 1986G 
and 2006ax at this time. 

4.4. Quantitative Measurements 

There have been several attempts to classify SNe la 
into distinct grou p s based on their spectra (e.g., 
Nugent et all Il995t iBenetti et all 120051 : [Branch et all 
2006|) . Most methods rely (at least in one dimen- 



sion] on the strength of Si II A5972 relative to that 
of Si II A6355. This is easily ni easure d by the pa- 
ramet er 7^(Si) (|Nugent et all 119951) . but [Branch et all 
(|2006f l advocates measuring t he equivalent widths of 
these lines. In addition to 7^(Si). [B"enetti et all ()2005l ) di- 
vided SNe la based on a combination of light-curve shape 
and the velocity gr adient of the Si II A6355 feature, v. 
iBranch et al.l (|2006l ) defines four subclasses, namely shal- 
low silicon (SS), core normal (CN), b road line (BL), and 
cool (CL), while iBenetti et all (|2005[ ) defines three sub- 
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Fig. 7. — Silicon ratio at maximum brightness (top; 7^(Si)max) 
and velocity gradient (bot tom; v(Si)) vs. A mi<j(B) for a sample 
of SNe la (taken from Be netti et al.l 120051) a nd SN 2006bt The 
LVG, HVG, and FAINT groups (as defined bv Benetti erani2005l ) 
are represented by squares, diamonds, and triangles, respectively. 
SN 2006bt is shown by the filled circle. Based on the velocity 
gradient, SN 2006bt is a member of the LVG group of SNe la 
despite having a spectrum similar to those of the FAINT group. 
SN 2006bt is a clear outlier in the relationship between the light- 
curve shape and the ratio of the silicon features. 



classes called low- velocity gradient (LVG), high- velocity 
gradient (HVG), and faint (FAINT). The classification 
systems have much overlap, with SS and CN objects 
mostly being in the LVG, BL objects mostly being in 
the HVG, and CL mapping directly to FAINT. 

At maximum light, SN 2GG6bt has W(5750) = 
29 A, 1^(6100) = 120 A, 7^(Si) = 0.44, and v = 
51kms~^ day^^. These values place SN 2006bt within 
the LVG group (although toward the high end of the 
range of v; see Figure [7|) and between the CL and BL 
groups, with its closes t neig hbors (using the updated 
values of iBranch erali 120091) being the BL SN 1999cc 
(VF(5750) = 26 A; 1^(6100) = 121 A) followed by CL 
SN 1989B (1^(5750) = 25 A; 11^(6100) = 126 A). De- 
spite being closest to a BL object in the equivalent width 
space, the shape of the Si II absorption and the border- 
line values should place SN 2006bt in the CL group. 

The true peculiar nature of SN 2006bt is apparent 
when comparing a temperature-dependent spectral pa- 
rameter such as 7?.(Si) with a light-curve parameter such 
as Ami^{B). In Figure [71 we show these values for 



the sample of [Benetti et a^ (|2005D (open symbols) and 
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Fig. 8. — Spectra of SN 2006bt (black curves) near the suspected 
C II AA6580, 7234 features. Spectra of SNe 2 006D (t = -5 days ; 
IThomas et al.l2007l) and 2006gz {t = -14 days: IHicken et al.ll2007l) 
are also shown for comparison (blue and red curves, respectively). 
Wavelengths near the potential C II features (corresponding to a 
blueshifted velocity of ~1350— 7850 kms~^) have been shaded grey. 
The C II features of SNe 2006D and 2006gz have blueshifted veloc- 
ities of 12,000 and 15,500 kms~^, respectively. [See the electronic 
edition of the Journal for a color version of this figure.] 



for SN 2006bt (filled circle). In this parameter space, 
SN 2006bt is a large outlier, indicating that it has a rel- 
atively cool photosphere and a slow decline rate. 

4.5. Possible Detection of Carbon 

In our earliest spectra, there is a small feature red- 
ward of Si II A6355 (see Figure [5]). In our first spectrum 
(when the feature is strongest and the signal-to-noise ra- 
tio of the spectrum is highest), its minimum is at 6465 A. 
This feature is present in our second and third spec- 
tra, although with less significance. These three spec- 
tra come from different instruments and telescopes, in- 
dicating that it is not an artifact of data reduction. The 
feature has mostly disappeared in our fourth spectrum 
{t = 0.0 days). 

An intriguing possibility is that this feature is from 
C II A6580 absorption. If this is the case, the 
minimum of the absorption would be blueshifted by 
5200 kms~^, which is much lower than that of other 
features and may indicate a misidentification. This fea- 
ture has bee n detected in a handful of objects (e.g., 
iBranch et aT 2003) with the be st cases being SNe 2 006D 
(IThomas et al.,.2007.). 200 6gz jHicken et alJ[2007h . and 



2009dc (jYamanaka et al.l l2009( ). but all detections so 
far have been at much larger velocities and typically 
at much earlier phases. For instance, the minimum of 
the C II A6580 feature in SN 2006gz was blueshifted by 
^15,500 kms^^ a.t t — —14 days (which was a larger 
velocity than the Si II feature at that phase) and had 



completely disappeared by i = — 12 days. 

C II has four relatively strong features in the optical: 
AA4267, 4745, 6580, and 7234. Unfortunately, C II A4267 
would be on the edge of a strong Mg II feature and 
C II A4745 is on top of a Ti II feature (although it may 
stiU be detectable with SYNOW; see Section [f6|. Only 
the two reddest lines are in relatively uncontaminated 
regions of the spectrum. As seen in Figure [H there is 
a very shallow feature in our earliest spectrum at the 
wavelength that corresponds to C II A7234 blueshifted 
by ~5000 kms"^ 

Observations have shown that high-velocity features 
are common m SNe la (iMazzaU et al.l I2005L e.g.,). 
The data can be explained by circumstellar interaction, 
dense blobs of material "detached" from the rest of 
the ejecta, or a diffe r ing covering fraction of the ejecta 
(iGerardv et al.ll2004 [Mazzah et al.ll2005t iTanaka et all 
[200^. If there is a single carbon blob ejected at an angle 
away from our line of sight, the absorption velocity seen 
in the spectra would correspond to the component of the 
velocity vector along our line of sight. For an angle of 45° 
(65°) relative to our line of sight, we infer a true velocity 
of 7450 kms~^ (12,300 kms~^) in such a scenario. This 
potential asymmetry may help explain other aspects of 
this SN (see Section [OJ. Alternatively, there may be 
a significant amount of carbon-rich material close to the 
SN. For this scenario, viewing-angle effects necessitate 
that our line of sight be in the plane of the circumstellar 
accretion disk. 

Besides the possible detection of carbon lines, 
SN 2006bt shares other characteristics with SNe 2006D 
and 2006gz. SN 2006D had a relatively narrow light 
curve similar to that of S N 1986G, but not as nar- 
row as that of SN 1991b g (IThomas et al.ll2007D . The 
t = —5 day spectrum from Thom as et al.l ( 20071 ) is repro- 
duced in Figure (5] The early-time spectra of SNe 2006D 
and 2006bt are relatively similar, both having relatively 
large 7^(Si). SN 2006bt has slightly larger velocities for 
the Si and Ca features and stronger Ti II features than 
SN 2006D. 

SN 2006gz had very broad light curves (Ami5(_B) = 
0.69 ± 0.04), no shoulde r in r, and no second maximum 
in i (IHicken et al. l2007i ). The colors were also very red, 
but IHicken et al.l (|2007f ) attributed this to a large amount 
of host-galaxy extinction {E{B — V)host = 0.18 mag). 
SN 2006gz had spectra indicative of a very hot photo- 
sphere unlike SN 2006bt. 

4.6. SYNOW Model Fits 

To investigate the details of our SN spectra, we 
use the supernova spectrum-synthesis code SYNOW 
(|Fisher et al.lll997D . Afihough SYNOW has a simple, 
parametric approach to creating synthetic spectra, it is 
a powerful tool to aid line identifications which in turn 
provide insights into the spectral formation of the ob- 
jects. To generate a synthetic spectrum, one inputs a 
blackbody temperature (Tbb), a photospheric velocity 
(uph), and for each involved ion, an optical depth at a 
reference line, an excitation temperature (Tcxc), the max- 
imum velocity of the opacity distribution (wmax), and a 
velocity scale (vg). This last variable assumes that the 
optical depth declines exponentially for velocities above 
Wph with an e- folding scale of Ve- The strengths of the 
lines for each ion are determined by oscillator strengths 
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Fig. 9.— Our -3.8 day optical spectrum of SN 2006bt (black) 
and our best-fit SYNOW synthetic spectrum with and without 
C II (blue dotted and red dashed lines, respectively). The insets 
show the regions near C II AA4745, 6580, 7234. The fit is improved 
slightly by including C II. [See the electronic edition of the Journal 
for a color version of this figure.] 



and the approximation of a Boltzmann distribution of 
the lower-level populations with a temperature of Toxc- 

In Figure [HI we present our —3.8 day spectrum of 
SN 2006bt with a synthetic spectrum generated from 
SYNOW. This fit is very simil a r to t he fit performed 
for SN 1986G by iBranch et all ()2006D and the optical 
depths for the fit are presented in Table [51 The main 
differences are that we include O I (it was excluded from 
the reference fit because of the short wavelength range of 
the SN 1986G spectrum), the excitation temperature is 
8000 K instead of 7000 K, all ions other than Mg I and 
Mg II have lower optical depths, we do not include high- 
velocity Ca II and Fe II, and we do not impose a maxi- 
mum velocity for Si II. The model includes low-ionization 
ions and excludes high-ionization ions, confirming the 
low-temperature nature of SN 2006bt. No ions at high 
velocity are included, further limiting the possibility of 
high-velocity ejecta. There are some differences between 
the SN 2006bt and SYNOW spectra. It is possible that 
these differences could be reduced by examining more 
of the parameter space of SYNOW input parameters. 
These differences do not affect our conclusions below. 

A slightly higher excitation temperature and lower op- 
tical depths in SN 2006bt compared to SN 1986G is con- 
sistent with our relatively crude assessment of the spec- 
trum in Section 14.21 where we considered SN 2006bt to 
be slightly hotter (and more "normal") than SN 1986G. 
The model for SN 1986G has Vn,^^ = 15,000 kms^^, 
while the model for SN 2006bt imposes no limit. Where 
nuclear burning may only extend to v = 15,000 kms~^ 
in SN 1986G, it appears that it extends to all velocities 
in SN 2006bt. 

Since we may have identified C II in the pre-maximum 
spectra of SN 2006bt (see Section|4?5|), we created models 
with and without C II. In an attempt to match the pos- 
sible low-velocity features, we included C II with Tcxc = 
12,000 K, Wmax = 15,000 kms-\ and = 200 kms'^. 
The inclusion of C II with these parameters improves the 
model fit. In addition to better fits near C II AA6580, 



7234, the model is significantly better near C II A4745. 
This shows that it is possible to detect lines in c onfused 
regions such as near 4600 A. (|Branch et al.ll2006[ ) fit this 
region with high-velocity Fe II, and again, we caution 
that the identification of C II is not definitive. 

5. SUPERNOVA ENVIRONMENT 

SN 2006bt occurred 33.7 kpc (7.8i?eff, where i?cff = 
6.4" is defined by the Petrosian radius where 50% 
of the fiux is contained within a circle of ra- 
dius -Roff) from the nucleus of its host galaxy, 
CGCG 108-013, an SO/a galaxy (see Figure 
SN 2006bt lies between its host and a nearby galaxy, 
2MASX J15562803-I-2002482 (SN 2006bt is 23.7" east 
of the nucleus of 2MASX J15562803-f 2002482; ahhough 
SN 2006bt lies between these two galaxies, we are con- 
fident that CGCG 108-013 is the host galaxy; see Sec- 
tion l4.ip . This galaxy has z = 0.0463, corresponding to a 
recession velocity t hat is 4250 kms~^ grea ter than that 
of CGCG 108-013 (jAbazaiian et all 120091). Examining 
SPSS images of the SN environment ( Abazaiian et al.l 



I2009D , there does not appear to be any tidal streams be 
tween the two galaxies. Although both objects are likely 
members of the Hercules supercluster, it is unlikely that 
they have interacted. There is another faint source 4.3" 
to the northwest (corresponding to 2.9 kpc at the red- 
shift of CGCG 108-013) that has colors consistent with 
a background gala xy (with SDSS photo-z of 0.36 ± 0.17; 
iCsabai et ai1[2003l) or a globular cluster. If it were a glob- 
ular cluster at the redshift of CGCG 108-013, it would 
have Mg — —12.1 mag, which is ~2 mag brighter than 
Gl in M31. It is therefore more likely to be a background 
galaxy. 

A sp ectrum of the nucleu s of CGCG 108-013 from 
SDSS ijAbazaiian et al.l l2009 h shows no emission lines, 
suggesting little or no recent star formati on. Using the 
single -stellar populations (SSP) models of lJimenez et al.l 
(|2004[ ). the spectrum of CGCG 108-013 is best fit by a 
population with a power-law initial mass function, an age 
of 14 Gyr, and a metallicity of 0.6 solar (see Figure fTTj) . It 
is worth noting that Jimenez et al. (2004) has no model 
with an age > 14 Gyr. Although such ages are limited 
by the age of the Universe, if the model ages are system- 
atically larger than true ages, we may not be probing 
the full range of realistic models. Despite the excellent 
fit with an old SSP model, we cannot rule out additional 
(low-level) star formation over the past several Gyr. 

Given the position of the SN, it is likely that the pro- 
genitor of SN 2006bt was a member of the halo popu- 
lation. If this were the case, then the progenitor would 
likely have formed before the majority of the stars in the 
bulge of the galaxy and with a lower metallicity. In this 
scenario, it is quite likely that the progenitor was of low 
metallicity and had an age of > 10 Gyr. 

Alternatively, the progenitor may have formed in a 
galactic disk or the bulge and was ejected outward from 
an interaction with a supermassive black hole at the nu- 
cleus of the galaxy. To reach 33.7 kpc, the progenitor 
would have to travel at 3300, 330, 33, and 6.6 kms"^ for 
10 Myr, 100 Myr, 1 Gyr, and 5 Gyr, respectively. From 
its position alone, it is very unlikely that the progeni- 
tor of SN 2006bt was a massive star (M > IOMq) that 
formed near the galactic center and was ejected into the 
halo. It is possible that the progenitor was a low-mass 
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Fig. 10. — {I^eft) SDSS gri false-color image of the position of SN 2006bt and its surrounding environment. The image is 4' x4' in extent. 
The yellow box is 1' on a side and centered on the SN position. The SN has a projected galactocentric distance of 33.7 kpc relative to 
its host galaxy, CGCG 108-013, the large galaxy on the northeast edge of the yellow box. The nearby galaxy to the west has a recession 
velocity 4250 kms~^ greater than that of the nominal host galaxy and does not appear to be interacting with the host galaxy, (right) 
Zoomed-in image of the SN position corresponding to the yellow box on the left (1' X 1'). The SN position is marked by a circle of 1" radius. 
There is a faint source 4.3" (corresponding to 2.9 kpc at the redshift of CGCG 108-013) to the northwest that is likely to be a background 
galaxy (see text for details). [See the electronic edition of the Journal for a color version of this figure.] 
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Fig. 11.— Spectrum of CGCG 108-013 (black) and the best-fit 
SSP model (t = 14 Gyr; Z = 0.012 = 0.6Zq; blue). 



star ejected in this manner, which would probably make 
it younger and have higher metallicity than other stars in 
the halo. Regardless, it is very likely that the progenitor 
of SN 2006bt was a white dwarf and SN 2006bt was not 



the result of some exotic core-collapse event. 

6. DISCUSSION 
6.1. Models 

It has long been believed that the vast majority of 
SN la light curves can be described by a single param- 
eter which directly correlates with the amount of ^^Ni 
produced in the explosion (with more ^^Ni correspond- 
ing to more lumino us eve nts having broader light curves). 
iKasen fc WooslevI (j2007[ ) explain this further by noting 
that objects with larger Mni are hotter, making the re- 
combination from Fc III to Fe II occur at a later time, 
an d thu s have a broader S-band light curve. 

Kasen (2006) found that the depth of the trough, the 
height of the second maximum, and the delay between 
the first and second maxima in the / band all increase 
with s^Ni mass. This finding would imply that SN 2006bt 
had a small ^^Ni mass (model /-band light curves with 
0.2 < A/ni < O.SMq show a similar behavior to the i- 
band light curve of SN 2006bt), which is in direct confiict 
with the decline rate of all bands. 

The prominence of the second maximum is more gen- 
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erally related to the mass of all Fe-group elements in 
the e.iecta , not just ^^ Ni and the isotopes into which 
it decays (jKasenI [ 20 06?) . The second maximum can be 
reduced by decreasing the progenitor metallicity, which 
decreases the amount of stable Fe-group elements ( "sta - 
ble Fe" ) produced in the explosion (T immes et al.l[2003[ ) ; 
however, the effect is at most ~10%. Since stable Fe is 
produced via electron capture, it is possible that the ex- 
plosion simply did not have the physical conditions nec- 
essary to form a significant amount of stable Fe through 
electron capture. Finally, one can smear out the sec- 
ond maximum in / b y mixing th e stable Fe and ^^Ni 
throughout the ejecta (|Kasenll2006f) . This will also cause 
the maximum of the / band to occur later, which may 
be the case for SN 2006bt. 

SN 2006bt may have been a very asymmetric explosion. 
Although we see no indication of high-velocity ejecta in 
any of the spectral features, there is the possibility that 
a carbon blob was ejected ~65° from our line of sight. 
SN 1999by, a low-luminosity event, showed a very polar- 
ized spectrum corresponding to an explosio n with an as- 
pheri city of 20% viewed along the equator (jHowell et al.l 
l2001t ). Asymmetric models show that peak luminosity 
depends on viewing angle; however, these models do not 
exhibit significant difference s in ionization w ith viewing 
angle (Kascn fc Plcwa 2007; Isim et al.ll2007l ). 

It is also possible that SN 2006bt rose very quickly 
to maximum, cau sing a high lu minosity with a smaller 
amount of ^^Ni (jArnettI fl98l . There are poor con- 
straints on an explosion date from pre-explosion images 
([Lee & Li 2006). Our photometry starts only a few days 
before B maximum, but it does not appear to be faster 
than other objects with its decline rate. Although a fast 
rise is still a possibility, it appears to be an unlikely so- 
lution. 

6.2. Implications for Cosmology 

Given its relatively normal luminosity, SN 2006bt- 
like objects would still be detected in high-redshift 
SN surveys, unlike other peculiar objec ts such as 
SNe 2002cx (iLi et al.l l2 003) and 2008ha (|Folev et all 
l2009aHValenti et al.ll2009| ). Furthermore, given its B-^ 
color curve, a SN 2006bt-like object with some host- 
galaxy extinction will simply look like a normal object 
with that extinction plus an additional Ay « 0.5 mag 
of extinction. At some point too much extinction would 
lead to the exclusion of these objects from a SN sam- 
ple; most cosmological analyses avoid SNe la having ex- 
tremely red colors. 

Additionally, at high redshift, the rest-frame r and i 
band will be redshifted out of the optical, requiring ob- 
servers to rely on rest- frame B — V alone to determine the 
quality of the light-curve fits and measure the extinction. 
At low redshift, one could search for objects with large 
Ar7ii5(i?) but no secondary maximum in i and exclude 
any such objects; however, at high redshift this is not 
feasible without observing in the near infrared {JWST 
and, depending on the ultimately chosen design, JDEM 
might be able to measure the rest-frame i band out to 
z sa 1). For an object at z > 0.3 (where i is redshifted out 
of the optical) , one would have to rely on a spectrum (in 
combination with a light curve) to identify SN 2006bt- 
like objects. However, with the Si II lines redshifting out 
of the optical window at z w 0.5 and the generally low 



signal-to-noise ratio spectra of high-redshift SNe la, it is 
difficult to detect spec troscopically pecul iar objects from 
high-redshift surveys (iFolev et al.ll2009bD . 

If SN 2006bt does come from a very old progenitor 
and the age of the progenitor is the main reason for its 
peculiarity, then SN 2006bt-like objects may not exist 
at high redshift. At z = 0.5 and 1, the Universe had 
an age of ^8.6 and 5.9 Gyr, respectively. If the de- 
lay time between formation and explosion is longer than 
those time s, SN 2006b t -like ob jects will not exist at those 
redshifts. iFolev et all (|2009bD found no SN 1991bg-hke 
objects in a sample of 118 high-redshift SNe la, indicat- 
ing that SN 2006bt-like objects are probably not a large 
fraction of the early-universe SN la population. 

We have developed a Monte Carlo simulation to as- 
sess the impact of contamination of a population of 
SN 2006bt-like objects in a SN la cosmological sample. 
For normal SNe la, we generate 1000 light curves over a 
range of light-curv e shape paramet ers from the MLGS2k2 
"early" templates (iJha et al.ll2007 ^ over < z < 0.8 (500 
for z < 0.1 and 500 for z > 0.1), with distance modulus, 
/itruc, calculated in a nominal flat AGDM {Um = 0.3, 
h = 0.65) cosmolo gy. We parameter ize the input extinc- 
tion by a i?y = 3.1 IO'DonneTl (I1994D law, with Ay drawn 
from the glos distribution fe.g.. lIIatano "eTaLlfrggl . The 
light-curve shape parameter A is drawn randomly be- 
tween —0.3 and 1.3. The sampling and signal-to-noise 
ratio is matched to the expected values for the Pan- 
STARRS Medium Deep Survey. 

Since SN 2006bt has very different light curves from 
those of other SNe la, the MLCS2k2 template hght 
curves had to be modified to create simulated light curves 
of SN 2006bt-hke objects. Specifically, the SN 2006bt- 
like light curves must be intrinsically underluminous and 
red compared to a nominal SN la with the same light- 
curve shape. To determine the appropriate modifica- 
tions, each individual light curve of SN 2006bt was fit 
separately. The UBVr light curves all had best fits 
of A « —0.15, while the i-band light curve was best 
fit with A = -0.38. The MLCS2k2 model hght curves 
were shifted so that the peak magnitudes were the same 
as those of SN 2006bt in each band. We created 100 
SN 2006bt-hke hght curves for < z < 0.8 (50 for z < 0.1 
and 50 for z > 0.1). 

All light curves were fit using MLCS2k2 (|Jha et all 
I2007D . which provided estimates of Ay, A, and fx. The 
fitting was performed with both the glos and flatnegav (a 
flat prior which includes negative values of Ay) priors. 
Fits of all light curves, including those for SN 2006bt-like 
objects, had a reduced ^ 2, indicating that all would 
nominally be included in a cosmological analysis. 

To demonstrate how MLCS2k2 interprets our 
SN 2006bt-like model light curves as a function of 
redshift, we produced noise-free light curves spanning 
the redshift interval < z < 0.8. The light curves 
were fit by MLGS2k2 with a flatnegav prior (in this 
signal-to-noise ratio regime, the prior is inconsequential 
and the uncertainty is dominated by the MLCS model 
uncertainty). The resulting distance modulus residuals 
and derived Ay are presented in Figure [121 The 
derived distance modulus residuals are close to zero 
near z = 0, but increase to ^^0.1 mag at z = 0.1. 
After that, the residuals generally linearly decrease with 
increasing redshift, deviating by ~0.5 mag at z = 0.8. 
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Fig. 12. — Distance modulus residuals {top; fiMhCS ~ Mtruc) 
and derived Ay (bottom) from MLCS2k2 as a function of redshift 
for the noise-free model light curves for SN 2006bt. The blue and 
red points represent objects with redshifts less than and greater 
than z = 0.1, respectively. The grey regions are the model un- 
certainties. There are several discontinuities at particular redshifts 
corresponding to either an observed-frame filter no longer matching 
a rest-frame filter (e.g., the observed V dropping out at z = 0.699; 
solid lines) or a transition from an observed-frame filter being K- 
corrected to a particular filter changing to a different one (e.g., 
observed i being K-corrected to rest-frame V switching to rest- 
frame S at 2 = 0.570; dashed lines). The black, blue, green, red, 
and orange lines represent the rest-frame UBVri filters undergoing 
transitions, respectively. [See the electronic edition of the Journal 
for a color version of this figure.] 



At redshifts where either the mapping of observed filters 
to rest-frame filters during the K-correction process 
changes or when a rest-frame filter is no longer used in 
K-corrections (the filter "drops out"), there are large 
discontinuities in the residuals. These discontinuities 
are partly the result of our model for SN 2006bt and 
the normal SN la template model and partly the result 
of having the observer-frame light curves derived from a 
single rest-frame filter. 

The derived value for Ay has a very similar (but in- 
verted) function with redshift. Since our model assumes 
Av = mag at all redshifts, the derived value corre- 
sponds to the difference between the true distance mod- 
ulus and the derived distance modulus, modulo the dif- 
ference in peak absolute magnitude for SN 2006bt and 
the canonical SN la with the same A. Although the dis- 



tance modulus residuals appear to be quite large, the 
distance modulus is compared to the input cosmological 
model. When fitting an observed sample, one can only 
look for outliers from the best fit, and the inclusion of 
SN 2006bt-like objects will influence such a fit. Further- 
more, a deviation of a few tenths of a magnitude from the 
residual of the best-fit cosmology is expected for several 
objects in a high-redshift SN sample. 

To examine the impact of SN 2006bt-like objects on 
cosmological parameter estimation, we defined 7 distinct 
samples. The first sample had no contamination from 
SN 2006bt-like objects. The remaining samples had 1%, 
5%, or 10% contamination at all redshifts, or only for 
z < 0.1 (consistent with a very long delay time for such 
objects). The distance moduli found using MLCS2k2 
were used to constrain the cosmological parameters flM 
and w with the assumption of flatness. We also combined 
the SN data wi th the prior from the b aryon acoustic oscil- 
lations (BAO: [Eisenstein et al.l[2005f ). further constrain- 
ing the cosmological parameters. Since we fit the light 
curves using two separate extinction priors, we have a 
total of 14 different estimates of ^Im and w with and 
without BAO constraints. The parameters are listed in 
Table O 

As SN 2006bt-like objects become a larger share of the 
SN sample, the estimate of w consistently increases (be- 
coming closer to 0). The value derived from the uncon- 



taminated sample is w 



-0.927t°:°?^ and 



-0.893 



+0.069 
-0.072 



for the glos and flatnegav priors, respectively, while the 
the 10% contaminated (at all redshifts) sample yields 



-0.813 



-1-0.065 
-0.066 



for the same priors. 



w = -0.725l°;°p and 
respectively. This trend is true if the SN 2006bt-like ob- 
jects are spread throughout redshift or restricted only to 
low redshifts, but which of these two samples has more 
bias depends on the prior. For the SN-only cosmological 
fits, the value of flM is typically biased to lower values, 
but only the fits with the glos prior are substantially bi- 
ased. The 1(7 contours for the uncontaminated and 10% 
contaminated samples are shown in Figure fT3l 

From our basic simulations, we see that SN 2006bt-like 
objects can have a large impact on derived cosmological 
parameters; however, if the fraction of SN 2006bt-like 
objects in the sample is low, the bias can be much smaller 
than other systematic biases (if 1% of the low-redshift 
sample in contaminated by SN 2006bt-like objects, w 
shifts by only 0.007). We have identified one object with 
these peculiar properties out of 185 objects in the CfA3 
sample, but we have not examined all objects with the 
same scrutiny. Within the CfA3 sample, 20% of SNe la 
have measured Ay larger than that of SN 2006bt, which 
is a very conservative upper bound to the fraction of the 
contamination of the low-redshift sample. 

7. CONCLUSIONS 

SN 2006bt is a unique object that defies the conven- 
tional wisdom of how SN la observables should corre- 
late. SN 2006bt has slowly declining light curves, but its 
spectra indicate a cool photosphere. The combination 
of the broad light curves and intrinsic red colors causes 
light-curve fitters to indicate that SN 2006bt is heavily 
reddened. The object is a slight outlier to the WLR. 
Although SN 2006bt is underluminous for its light-curve 
shape, light-curve fitters correct for this in a way that 
estimates a distance modulus that is very close to the 
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TABLE 3 

COSMOLOGICAL PARAMETER RECOVERY FROM SIMULATIONS 



Sample 


(SNe Only) 


(SNe + BAO) 


W 

(SNe Only) 


W 

(SNe + BAO) 




flatnegav Prior 


No contamination 

1% contamination 

5% contamination 

10% contamination 

1% contamination low-z 

5% contamination low-2 

10% contamination low-2 


n Q1O+0.097 
U.010_Q 203 
O99+0.097 
u.ozz_Q 206 
314+0106 
u.ol^_Q 202 

27fi+° "'^2l 

U.Z(D_Q j^g 

n 900+0.139 

U.^3«_Q igg 


292+°-°28 
u.zyz_Q gig 

r. 900 + 0.028 
U.Z»0_Q gig 

296+° °28 

U.ZMD_g gig 

U.OUZ_g gig 

900+0.028 
U.Z»0_g gig 

U.Z»^_g gig 

296+0 028 

U.ZMD_g gig 


-1 079+0-346 

l.UI 3_g 47g 

-1 076+0-347 

i.U/D_g 47g 

1 m 0+0.322 

l.U10_Q 4gg 

-1 014+0-334 

l.U11_g 4g7 

1 0-^7+0-326 

l.UOf _g 4g2 

-0 921+0-266 

U.»Zl_g 4yi 

-0 795+0-200 


n OQq + 0.069 

U-030_g 072 

n SQQ+0.069 
-U.»»3_g 071 

«c:7+0.067 

U.OO/_g ggg 

n SI -3+0.065 

U.010_g Qgg 

-0.8861H?? 
-0-8641HS 

n s'?9+0-064 

U.OOZ_g Qgg 


glos Prior 


No contamination 

1% contamination 

5% contamination 

10% contamination 

1% contamination low-z 

5% contamination low-2 

10% contamination low-z 


n 01 1+0.082 
U.011_Q ig2 

0.278l°:f8j 

0.206l°;l|| 
,„_+0.164 

0.iSU_o 120 
n 151+0.171 


fj 9OQ + O.O27 
U.28y_g gig 

295+0-029 

U.ZMO_g gig 

302+0 029 

U.OUZ_g gig 

n qnS+0 029 
U.3U«_o g2o 

295+0 029 

U.ZyO_Q gig 

302+0 029 

U.OUZ_g gig 

n 01 1 +0.029 

U.011_Q^g2g 


-1 094+0-324 
l.Ua4_Q 415 

-0 Q21+0-268 
"•"^^-0.476 
n fi7«+0-144 

— 61 Q+0.123 

U.Diy_g 335 

-0.852tO-230 
-0.633tO;^l? 

n K97+0.076 
U.OZ(_g 229 


-0 Q27+0-067 
^•az'_o.070 
n S61 +0.067 

-U.»bi_g ggg 

-0 774+0.060 
' '^-0.060 
„ 79P.+0.055 

U. (ZO_g 

S/18+0-065 

-0.848_g gg7 

-0 765+0-057 
U. IDO_g 

-0 678+0-051 


Note. — BAO alone results in Om = 0.293tQ 077 ^^'^ ™ 




Since the input 



cosmology of Q.m — 0.3 and w — —1 differs slightly from this best-fit result, there is some tension 
between the datasets and their combination will cause the combined contours to move toward 
the BAO-only results. 



-1.1 


10% z < 0. 1 

10%z<0.8 




-1.0 






-0.9 






S 

-0.8 






-0.7 






-0.6 







0.30 



0.30 



Fig. 13. — la confidence contours in the Qm^w plane for a 
combination of measurements of 1000 simulated SN la light curves 
with BAO measurements. The black solid line and dark grey region 
correspond to a sample of SNe la with no SN 2006bt-like objects, 
while the blue dotted line (and light grey region) and the red dashed 
line (and medium grey region) corresponds to samples with 10% 
contamination of SN 2006bt-like objects in the nearby (z < 0.1) 
and full sample, respectively. The green dot represents the input 
cosmology for the simulation {Qm = 0-3 and w = —1), and is 
offset from the most likely value of the uncontaminated sample 
because of the combination with the BAO measurements, which 
prefer Qm = 0.293 and ui = —0.920 (marked by the yellow star) 
before combination with SN data. The left and right panels are 
for light-curve fits using the flatnegav and glos extinction priors, 
respectively. [See the electronic edition of the Journal for a color 
version of this figure.] 



Hubble-flow value. 

We discovered the peculiar nature of SN 2006bt by 
noticing that it had a large d erived extinction value for 
its PGCD, but this was after |H09b| included SN 2006bt 
in a cosmological analysis. We may find similar objects 
by looking for comparable outliers, but it is also possible 
to find SN 2006bt-like objects by examining their light 
curves and spectra. SN 2006bt is a clear outlier in the re- 



lationship between Atois and TZ{Si). Unlike SNe la with 
similar decline raters, SN 2006bt lacks a prominent sec- 
ond maximum in the i band. Both of these observables 
could be used to find additional examples. 

There are spectroscopic features which are consistent 
with C II, but at a significantly lower velocity than other 
spectral features. If these features are caused by carbon, 
it can be explained by a carbon blob ejected ~65° from 
our line of sight. 

The SN position and the host-galaxy type and spec- 
trum indicate that SN 2006bt had a very old progen- 
itor. If a very old progenitor is necessary to create a 
SN 2006bt-like event, we may not expect to find any 
similar objects at high redshift. 

Contamination of SN 2006bt-like objects in samples 
of SNe la used for cosmological analyses is a potential 
source of a large bias in the determination of and 
w. If these objects are only ~1% of the lo w-redshift 
sample (similar to the sample used by iHOQbT ) , then the 
bias should be < 0.01 in w. However, if these objects 
are a larger percentage of the sample, this could be the 
dominant systematic uncertainty for SN cosmology. It 
seems unlikely that SN 2006bt-like objects are more than 
a few percent of the local sample, but it is much harder 
to characterize similar objects at high redshift (where 
the rest-frame / band is not observed, and low signal-to- 
noise ratio spectra may not cover the Si II A6355 feature) . 
Future studies should attempt to identify SN 2006bt-like 
objects and exclude them from cosmological analyses. 

Facilities: FLW0:1.5m(FAST), Keck:I(LRIS), 
Lick:Shane(Kast), Sloan Digital Sky Survey 
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